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Abstract 

We study a disk amplitude which has a comphcated heterogeneous matter configuration 
on the boundary in a system of the (3, 4) conformal matter coupled to two-dimensional grav- 
ity. It is analyzed using the two- matrix chain model in the large N limit. We show that 
the disk amplitude calculated by Schwinger-Dyson equations can completely be reproduced 
through purely geometrical consideration. From this result, we speculate that all heteroge- 
neous loop amplitudes can be derived from the geometrical consideration and the consistency 
among relevant amplitudes. 
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1 Introduction 



The {m, m + 1) unitary minimal conformal model coupled to two-dimensional gravity can be 
analyzed by the two-matrix chain model or the multi-matrix chain model |l|-||5| under appro- 
priate scaling limits. In two-dimensional gravity, loop amplitude is one of the most important 
physical quantities, which is the amplitude for a surface with loop-like boundaries. The (m, m+1) 
minimal model has m — 1 matter states [^], so in general, we have to deal with loops on which 
matter configurations are not restricted to homogeneous ones [^-[12|. 

In Refs. |11, 12|, some disk amplitudes whose boundaries have heterogeneous matter states 
are studied in the case of the tricritical Ising model (m = 4 case). It is found that the obtained 
amplitudes have simple geometrical interpretation. The heterogeneous loops make interesting 
splitting and sticking interaction regularly. Here a question arises. Can we decide any het- 
erogeneous loop amplitudes only from the geometrical consideration? In other words, does the 
geometrical consideration give enough information to decide any heterogeneous loop amplitudes? 

In order to address this question, in this paper, we study a disk amplitude with rather 
complicated boundary condition in the case of the Ising model (m = 3 case). We show that 
the amplitude calculated by Schwinger-Dyson equations can be reproduced through purely ge- 
ometrical consideration. We speculate that any heterogeneous loop amplitudes can be decided 
geometrically. 



2 Amplitude from Schwinger-Dyson equations 

In this section, we examine a disk amplitude with heterogeneous matter boundary condition 
using Schwinger-Dyson equations. Let us start with the action of the two-matrix model 

S{A,B) = ^ti{U{A) + U{B)-AB}. (2.1) 

Here A and B are N x N unitary matrix variables, and A is the bare cosmological constant. As 
a critical potential which realizes the Ising model coupled to two-dimensional gravity, we take 

[/((/<) = «(/.+ 1^2 + 2c/>3 , (a, A 7) = (3,-3,-1) . (2.2) 

This can be found using the orthogonal polynomial method |5|, pT|] . 
In this paper we will focus on the amplitude 

oo . 

WABAB{pi,qi,P2,q2) = E {tT{A^ A^B")) pi'''\^^-'p^^^^-\,'^-' (2.3) 

k,l,m,n=0 

and its continuum universal counterpart WABAsiCAi, Cbi, Ca2, (32,^) in the large limit. Here 
Pi and Qi are bare boundary cosmological constants, Ca^ and ^Bi are their renormalized counter- 
parts, and t is the renormalized cosmological constant. The insertion of the operator ti {A^B''A"^B^) 
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makes a loop-like boundary on which a part with spins up and a part with spins down appear by 
turns. In the large N limit, the corresponding continuum amplitude becomes the disk amplitude 
whose boundary consists of four arcs with distinct matter states. 
Let us start with the Schwinger-Dyson equation 

= J2 f dAdB-^ {tr [aH^B^A^B"^^ e-S{^,B)| ^ (2.4) 

a •' 

here the hermitian matrix A is expressed as A = J2a A°'t°' using the bases of hermitian matrices 
{f^}- Rewriting Eq. (2.4) explicitly, we obtain 

= + [B^A"'-^-^][A''A^B''] - a[A^B^A"'B''] 

i=0 i=0 

-P[A^+^B^A"'B''] - 7[^'=+25'A"S"] + , (2.5) 

here we use the notation [A'^B^---] = ^(tr(^^S' • • •)) • Multiplying Pi'''^qi^'^P2"^'^q2 '^'^ 
and summing up with respect to k, I, m and n indices, we obtain the equation in terms of 
resolvent amplitudes 

= (WAipi) -viPl) + qi) WABAB{pi,qi,P2,q2) + (WBAiqi,P2) - 1) WAAsiPl, P2, q2) 

+ {P + lPl)WBAB{qi,P2,q2) + iW^"^^ BABiqi,P2,q2) ■ (2.6) 

In Eq. (^) and in the following, we use the notations v{p) = a + Pp + jp"^ , 

OO A 

I^M.M....(n,r2, •••)= E ^(tr(MfM^••))r^'=-V2-'-^ ... , (2.7) 

k,l,-=0 

OO 4 

T^(*'V(*''V....(ri,r2, •■■)= E ^(tr(MMl^M'M^..))r^'-V2-'-l ••• , (2.8) 

k,l,---=0 

where M, M' and Mj denote the matrices ^ or i?. Using the relations WBA{qi,P2) = W^ab(P2; 9i) 
and WBAB{qi,P2,q2) = WAAB{q2,qi,P2), we observe that Wabab can be expressed in terms of 
Wa, Wab, Waab and W'^'^^bab- 

As for Pr(^) BAB 1 from the resolvent expression for the / = case of Eq. ( |2.5| ) and the relations 

W'^'^'>ABiP2,q2) = P2WAB{P2,q2) - ^^5(92) and Wa^^'' Ab{P1, P2, q2) = W ^^'^ b AB iP2 , q2 , Pi) , We 

obtain 

W^"^^ bab{P2, q2,Pi) = - {Wa{pi) + Wa{P2) - v{pi)) Waab{pi,P2, 92) 

-{fi + l{Pl+P2))WAB{p2,q2)+lWB{q2) ■ (2.9) 



Combining Eqs. (|2.6|) and (|2.9D , we find that Wabab can be expressed in terms of Wa, Wb, 
Wab and Waab- Furthermore using the relations 

TT. I N WAB{p\,q) -WABij>2,q) 

WaabKV\,P2, q) = , (2.10) 

pi -P2 
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Wa{p) - v{p) + q 

WB^^\q) = - (Wsiq) - viq)) Wsiq) - (/? + 7'?)A - -f[A] , (2.12) 

we find, in the end, that Wabab can be expressed in terms of Wa and Wb- (Here the last two 
equations are obtained in a similar way to Refs. 0, |^, 11, |l2|.) 

The homogeneous disk amplitudes Wa and Wb satisfy a third order equation. They are 
expanded as 

/„\ 4/3 

WA{p) = WB{p) = 3-2aC+i^-j w{C) + 0{a^/^) , (2.13) 
here a is the lattice spacing and 

w{c) = {c+^/e^t f\[c-^/e^tf (2.14) 

is the continuum homogeneous disk amplitude ||l^ under the renormalization A = 10 — a'^t and 
p = aC- This amplitude satisfies the third order equation 



w 



{Cf - ?>t^'^w{0 - 16C^ + - 2*2 = . (2.15) 



Combining all of the above equations, we obtain the continuum heterogeneous disk amplitude 
WABAB- From the explicit calculation, we find that terms in 0{aP^^) in Eq. (|2.13D do not affect 
the calculation of wabab-, and wabab is expressed only in terms of w{C,Ai) w{C,Bj)- From 
the a} order terms of 16 Wabab, we obtain the continuum disk amplitude wabab as follows: 



WABAB (Cai , Cbi , Ca2 ) Cb2 , 

w{C,A^)^w(C,A2) i w(C,A^'^ w{Ca2)' 



CAi - CA2 V (^^1 + CBi)(CAi + (,B-i) (CA2 + CBi)(CA2 + Cb2 



CAi - Ca2 V (^^1 + ''^i ) (^^1 + ''^a ) (Ca2 + Cbi ) (Ca2 + Cb2 ) 
w{C,B^)w{C,B2) ( w{C,Ar) W^Ca^) 



CAi - CA2 V (CAi + CBi ) (CAi + Cb2 ) (CA2 + CBi ) (CA2 + Cb2 ) 

CAi + CA2 + CSl + Cb2 



-3t^/3 (.(cJ + .(Ca2)) (c^CbJ(Ca2 + C.J(Ca. + Cb2)(Ca2 + CB2) 

+ (Ai ^Si, ^^2). (2.16) 

Here the last expression (Ai Bi, A2 <-> -B2) denotes the terms which are obtained by inter- 
changing Ai and Bi in all of the preceding terms. As the continuum part, we have taken terms 
which are non-analytic in all of the variables Ca^) Cbj and t. We should notice that the form 
of the amplitude ( 2.16| ) is not unique because of the identity equation ( p. 15 ). This uncertainty, 
however, does not arise in the geometrical consideration in the next section. 
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3 Amplitude from geometrical consideration 



In this section, we will show that the heterogeneous disk amplitude wabab obtained in section 
2 can be reproduced through purely geometrical consideration. 

At first let us briefly review the mechanism of the interaction of heterogeneous boundary in 
the case of wab IHIi which is the continuum counterpart of Wab- In this case, the boundary loop 
consists of two parts on which the spin state is up or down. According to boundary conformal 
field theory, discontinuity in the matter boundary conditions corresponds to the insertion of a 
boundary operator in Refs. [0, |9|. So in the case of wab, we consider that two spin operators 

(f)^ = a are situated at the points where spin state changes (see Fig. ||). The continuum 

amplitude wab is calculated to be ||9|, |l^ 

wab{Ca,Cb) = 7 — — ^ • (J-lj 

By performing the inverse Laplace transformation, we obtain the disk amplitude VV^_b(^A)^b) 
in terms of the lengths of the boundaries (see appendix A), 

WABi^AjB) = Ca-^Cb-^[wab{Ca,Cb)] 



Q{U - ^b)(W * W)(^A - ^b) + G{tB - iA)^ * >V)(^ 
+ / d£ W{iA - W{£b -i) - 3t^/^S{£A - £b) ■ (3.2) 







From Eq. (3.2) we find that the following geometrical interpretation is possible |11]. The orig- 
inal heterogeneous loop splits into two homogeneous loops. The first three terms in Eq. ( |3.2| ) 
correspond to the three diagrams in Fig. |^ respectively. This loop splitting phenomenon can be 
summarized in the following rules. 

1. Boundary A |^ sticks to boundary B so that a heterogeneous loop splits into homogeneous 
loops. The split disks are connected by the double line which is formed by sticking of 
boundaries A and B. 

2. The spin operator a is situated at each end of the double line. 

The spin operator a separates the boundary into two sides. One side forms homogeneous loop- 
like boundary and the other forms a stuck double line. We speculate that the spin operator a 
has the effect to make discontinuity in the geometrical state of boundary as well as in the matter 
state. 

Now let us derive wabab through a purely geometrical consideration. For simplicity, we 
will neglect the terms proportional to t^/^ in Eqs. ( 2.16| ) and (3.1), because they only make 



^ We refer to the part of boundary on which the spin state is up as boundary A etc, in the same way to p^. 
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the argument complicate. We can deal with these terms in a similar manner. The rules of the 
splitting phenomenon should equally be applied in this case. We will employ the rules and sum 
up all possible terms. Here the weights of the respective terms, however, are not fixed. The 
relative weights will be fixed by requiring the consistency with wab and the expected symmetries 
of the amplitude wabab- The resultant amplitude will completely coincide with that derived 
from the Schwinger-Dyson equations up to overall normalization. 

In the following we will explain how to derive the amplitude wabab in detail. In the case 
of Wabab, there are four spin operators a on the loop- like boundary. Let us focus on one of 
them (denoted by o in Fig. After loop splitting, the spin operator we focus is in touch with 
another part of the boundary. According to the point with which it is in touch, there are four 
possibilities, (a), (b), (c) and (d) in Fig. ^. As an example let us consider the case (a) concretely. 
The spin operator we focus is at one end of the double line (see Fig. ^). At the other end of the 
double line, one of the remaining spin operators will be situated. According to the spin operator 
situated there, the case (a) is classified into three sub cases, (a- 1), (a-2) and (a-3) in Fig. |. 
Then the amplitude corresponding to the case (a) is written as follows (see appendix A): 

g ^ ^AAB (CAi , CAz , CBi ) w{Cb.2 ) ^ WBBA (CBi , Cbz , CAi ) w{Cb2 ) 
CA2 + Cb2 CA2 + Cb2 

{Q3w(Cai) w(Ca2) +Q4w(Cai) w(Cbi) + a5w{CB2)w{C,A2) + a'iiW{CB2)w{C,Bi)}w{CB2) 

(CAi + Cb2 ) (CA2 + CBi ) (CA2 + Ci?2 ) 

(3.3) 



+ 



Here Oj represent unknown weight constants, which will be determined later, and waab and 
WBBA represent the continuum counterparts of the amplitudes Waab and Wbba respectively. 
Similarly the case (b) is classified into two sub cases, (b-1) and (b-2) (see Fig. ^). The corre- 
sponding amplitude is 

^ WAAB (CAi , CAz , Ci?i )w{C,B2) ^ WAAB (CAi , CAz , Ci?i ) ^^(CAi ) 
CAi + Cb2 CAi + Cb2 

where hi are unknown weight constants. The cases (c) and (d) are obtained by interchanging 
with Bi in the cases (b) and (a) respectively: 

Qc = Qb\A^B-, Qd = Qa\A^B- (3.5) 

In the cases (b) and (c), we have doubly counted several terms. The next amplitude corresponds 
to the doubly counted terms, 

r. _ u^^^m)w{C,B2){w{Ca-,) +w{Cb^)) ^^w{Cai)w{Cbi){w{Cai) +w{Cbi)) 



(CAi + Csi ) (CAi + Cb2 ) (Ca2 + Csi ) (CAi + Cbi ) (CAi + Cb2 ) (Ca2 + Cbi ) 

^ w{CBi)w{CB2){eiw{CAi) +e2w{CBi) ) ^ w(Cai)w(Ca2) (ei u'(Cbi) +e2w(CAi) 



(Cai + Csi ) (CAi + Cb2 ) (Ca2 + Csi ) (CAi + CBi ) (CAi + Cb2 ) (Ca2 + Cbi ) 

(3.6) 
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Here the first term comes from the double counting between the cases (b-1) and (c-1) (see Fig. ^). 
The second term comes from that between the cases (b-2) and (c-2) (see Fig. ^). The third term 
represents that between the cases (b-1) and (c-2) (see Fig. ^). The last term represents that 
between the cases (b-2) and (c-1) (see Fig. |lO|). 

The amplitude we aim at is obtained by summing up all the cases 

G = Ga + Qb + Qc + Qd + Ge ■ (3.7) 

By using the relations 

WAABiCl,C2,0 =WBBAiCl,C2,0 = 7 7 , (3.8) 

WABiC,^) - T—z , (3.9) 

c + ? 

we can express wabab in terms of the homogeneous disk amplitudes. It contains 10 unknown 
weight constants. They can be decided by requiring symmetries of the amplitude. Eq. ( |3.7| ) 
already has the symmetry under interchanging Ai with Bi. We further require the symmetry 
under (Ai <-> A2) and that under (Bi ^ B2) ■ The unknown constants are determined up 
to overall normalization. The weights (01,025041061^2) have the same non-zero value and the 
remaining constants equal to zero. 

In this way, wabab is determined up to the normalization constant. When we choose the 
normalization constant properly (oi = 1) , the obtained amplitude indeed coincides with that 
obtained from the Schwinger-Dyson equations in section 2 R. 



4 Discussion 

In the previous section we derived the heterogeneous disk amplitude wabab only through 
geometrical consideration. We required the symmetries of the amplitude there, however, we 
consider that the rules of splitting and the consistency with wab were essential. 

Let us discuss to what extent we can generalize the result we obtained. According to the 

rules of loop splitting, we can expect that wab consists of three terms^see Fig. §). The rules, 

however, do not give the information on the relative weight of each term. That is, the weight 

of the last term of the numerator of Eq. ( |3.9[ ) relative to the first two terms is unknown. When 

we treat this weight as an additional unknown constant, we can also determine it in the same 

way to section 3. This shows that we can also derive wab from the rules of loop splitting and 

the consistency between wab and wabab- From this fact, we speculate that any heterogeneous 

^ In a similar manner we can also reproduce the terms proportional to t'^^'^ in Eq. ( |2.16| ). 
^ Here we neglect the term proportional to t'^^^ . 
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loop amplitudes can be determined from the loop-splitting rules and the consistency among the 
relevant amplitudes. 

In boundary conformal field theory, the boundary operator a plays the role to make discon- 
tinuity in matter state on the boundary. As we saw in section 3, when it is dressed with gravity, 
it also makes a discontinuity in the geometrical state of the boundary. Here let us comment on 
the relation with the boundary operators ^ which appear in the scaling operators in the 
matrix models. In Ref. [16| it is argued that the scaling operator (T„(m+i) = has the effect to 



split a loop into k loops in the case of (m, m + 1) model, where k < n. When two a approach 
each other, one can consider that they change the spin state on the boundary locally, and they 
have the effect to split a loop into two loops |9|, |l2|. The boundary operator B2 can be identified 
with the operator obtained in the limit where two a approach each other. We expect that i?„ 
can be constructed by making a linear combination of the operators obtained when 2{k — 1) spin 
operators a get together, where 1 < k < n. 
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Appendix A Inverse Laplace transformation 

We summarize here some useful relations on inverse Laplace transformation. The Laplace 
transformation of a function f{ii) is defined by the relation 

Ci[fm = F{Q) = / dii e-'^^^ . (A.l) 







The following three relations are useful in calculating the inverse Laplace transformation, 

1 



r-lr-l 



'-1 '-2 



C1 + C2 
^(Ci) 



5(4 - 4) , (A.2) 
- ^2) /(4 - ^2) , (A.3) 



.C1 + C2 

A"' [^(C.) cm = t M' /(/) 5(4 - /) = /(4) *i g{^^) . (A.4) 

Here -F(Ci) and G(Cj) represent the Laplace transformed functions of /(4) and g{ij) respectively, 
and the symbol *j in the third equation denotes the convolution with respect to 4- All of the 
inverse Laplace transformations of the amplitudes in section 3 can be performed by combining 
these relation^. 



Eq. ( |A.3[ ) can be derived combining Eqs. (A..2) and (A.4) 
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As examples, we show the inverse Laplace transformations of some terms in Eq. ( |3.3| ) in the 
following: 

" WAAB (CAi , Ca2 , Cbi ) w{Cb2 ] 



n-\ p-\ p-\ 

'-Ax^A^'^Bx ^B2 



CA2 + Cb2 

WA2-^B2)y^AAB{Ux.^A2- 

WAAB(.iA„iA2 - i'jB,) W{£b2 - i') , 



„^bJ] *B2 W(^B2) 



(A.5) 



r-l r-l r-l r-l 
'-Ax'^A'i '-Bx'-Bi 

= Wa, - 

*B,{[Oi 



w{Cax) w{Ca2) w{Cb2) 

(CAi + Cb2)(CA2 + CSi)(CA2 +CB2). 

iB,)mA, - iB,)] 

U2 - ^bJW(£a2 - iB,)] *A, Wb, - iA,)W{lB, - £a,)]} . (A.6) 



From Eq. ( |A.5| ), one finds that this term corresponds to the first diagram in Fig. ^. Similarly 
the amplitude (|A.6| ) can be represented by Fig. 11. 
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A2 

Figure 3: There are four cases 
depending on the behavior of the 
spin operator we focus. 




Figure 4: This shows the interme- 
diate situation corresponding to 
the case (a). Observe that the 
spin operator we focus is situated 
at the right end of the double Une. 




(a-3) 

Figure 5: The case (a) is furthermore classified into the three sub cases. 
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Figure 7: The doubly 
counted diagram between 
(b-1) and (c-1) 
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Figure 9: The doubly 
counted diagram between 
(b-1) and (c-2) 



Figure 8: The doubly 
counted diagram between 
(b-2) and (c-2) 



B1 or A1 




Figure 10: The doubly 
counted diagram between 
(b-2) and (c-1) 
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